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Abstract 

Based on an effective Lagrangian which combines chiral SU(3) dynamics with 
vector meson dominance, we have developed a model for the forward vector meson- 
nucleon scattering amplitudes. We use this as an input to calculate the low energy 
part of the current-current correlation function in nuclear matter. Its spectrum en- 
ters directly in the "left hand side" of QCD sum rules. For the isovector channel we 
find a significant enhancement of the in-medium spectral density below the p reso- 
nance while the p meson mass itself changes only slightly. The situation is different 
in the isoscalar channel, where the mass and peak position of the a; meson move 
downward while its width increases less drastically than in the p meson case. For the 
(j) meson we find almost no mass shift; the width of the peak broadens moderately. 
We observe a remarkable degree of consistency with the operator product expansion 
of QCD sum rules in all three channels. We point out, however, that these results 
cannot simply be interpreted, as commonly done, in terms of a universal rescaling 
of vector meson masses in matter. 



Introduction 



At present there is a lively discussion about the behaviour of vector mesons in dense 
and hot hadronic matter. According to Brown and Rho vector meson masses my (in 
fact the masses of most hadrons with the exception of the pion) should drop in matter 
according to a simple (BR) scaling law: my/mv = f^/ fn, where /^r is the pion decay 
constant and the asterix refers to in-medium quantities. A decreasing /* (more precisely; 
the /* related to the time component of the axial current) in matter indicates a tendency 
toward chiral restoration (the transition from the spontaneously broken Nambu-Goldstone 
to the unbroken Wigner-Weyl realization of chiral symmetry). Dropping vector meson 
masses would then be a signal for this transition, if BR-scaling holds. Other scaling laws 
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0] using a bag model find a mass reduction of two tliird times tliat of the nucleon mass 
which leads to a similar shift as BR-scaling. 

Several analyses of the p and uj meson masses in matter using QCD sum rules 0, |^ 
seem to confirm these expectations. On the other hand, model calculations of the density 
dependent two-pion self-energy of the rho meson in nuclear matter 0, ^ suggest only 
marginal changes of the in-medium rho meson mass, but a strongly increased decay width 
instead. Additional BR scaling seems to be needed in order to match the hadronic models 
with the QCD sum rule analysis 0. Phenomenologically, a dropping p meson mass seems 
to help 1^, 3 understanding the enhanced dilepton yields seen at masses below the rho 
meson resonance in the CERES [[l^] and HELIOS |lT| experiments at CERN. 

Few attempts have been made so far to perform a realistic calculation of the omega 
meson self-energy in matter. The influence of the density dependent effective nucleon 
mass in the particle-hole polarization coupled to the uj meson has been studied in ref . Jl^ . 
However, leading in-medium effects from the cu-decay into the three-pion continuum in 
the presence of nucleons have not yet been investigated. One of the aims of this paper is 
to take such effects into account. 

The meson as the third neutral member of the SU(3) vector meson octet is special 
in the sense that it is almost a pure strange quark- ant iquark system. The large strange 
quark mass protects it from effects related to chiral restoration. Therefore the meson 
mass does not seem to be rescaled in matter as strongly as the p oi uj meson when analyzed 
using in-medium QCD sum rules. 

The QCD sum rule approach, at any practical level so far, has used only a caricature of 
the true spectral functions, namely a parameterization in terms of a delta function at the 
vector meson pole accompanied by a step function continuum at higher energies. While 
such a parameterization turns out to be valid in vacuum (as we shall confirm), this is shown 
not to be justified for the in-medium spectra. In particular, the parameterization of the 
rho meson in matter as an isolated quasiparticle pole becomes unreliable with increasing 
baryon density, given the strong increase of its partial decay width into pN — > vrA^, ttttA^ 
etc. channels. 

The purpose of this paper is to present a systematic calculation of current-current 
correlation functions in baryonic matter, based on an effective Lagrangian which combines 
chiral symmetry with vector meson dominance of the electromagnetic current. In the 
isovector channel we find a significant enhancement of the low mass part of the spectrum 
with increasing baryon density. On the other hand the p meson mass changes at a much 
slower rate with density than naively expected. In the isoscalar channel the situation 
turns out to be different. The position of the uj meson pole moves downward with rising 
density. The uj width does not increase as strongly as for the isovector channel. To leading 
order in the density, these features can be translated into statements about the (complex) 
effective scattering lengths of p and uj interacting with a nucleon. These will be discussed 
in detail. The meson-nucleon scattering amplitude is much smaller than the pN and 
ujN amplitudes, with a slight increase of the width but almost no shift of the resonance. 
We show that the vector meson-nucleon scattering amplitudes derived from our hadronic 
effective Lagrangian lead to spectral densities which are remarkably consistent with the 
standard QCD sum rule analysis, without the need for additional BR scaling. 

This paper is organized as follows. In section we review properties of the current 
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correlation functions in vacuum and give a brief account of the QCD sum rule approach. 
The current-current correlation function in baryonic matter will be developed in section 
^. To leading order in the density this requires a detailed calculation of the vector meson- 
nucleon scattering amplitudes which will be described in section ^. The in- medium spectra 
and their comparison with QCD sum rules is presented in section and the results will 
be discussed in section |[ 



2 Current-current correlation functions in vacuum 
2.1 Basics 

This is an introductory reminder of properties of the current-current (CC) correlation 
function in vacuum (for a review see e.g. ref. |]TB|). The vector mesons are seen as 



resonances in the current-current (CC) correlation function, 



where T denotes the time-ordered product and is the electromagnetic current which 
can be decomposed as 

J, = jl!'^+j!r^+jlf^ (2) 

into vector currents specified by their quark content: 

= ^{u-ff,u- d-f^d), (3) 
= l{u^,u + d^,d), (4) 
Jlf^ = -^is^.s). (5) 

Current conservation implies the following transverse tensor structure: 

n,.(g) = (^g,u - ^) n(g2), (6) 
with the scalar CC correlation function 

n(g^) = Ig^'^n.Aq). (T) 

Its imaginary part is proportional to the cross section for e^e~ hadrons: 

„ , , (7(e'^e~ — i> hadrons) 127r^ , 

Ris) = = imn s , 8 

cr(e+e fi^fi ) s 

where y/s is the total cm. energy of the lepton pair and a{e~^e~ /^^A^ ) = 47ra^/3s 
with a = e'^/Air = 1/137. The vector mesons can be distinguished by looking at different 
hadronic channels carrying the respective flavour (isospin) quantum numbers (see eqs.(B]- 
^). The corresponding correlation functions are denoted by 11^^^''^ 11'^^"^'' and 11'^'^. 
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G-parity demands that the isovector current involving the p meson can only couple to 
even numbers of pions (see the data in fig. la). Similarly, the annihilation of e^e~ into 
odd numbers of pions determines the isoscalar current involving the uj meson, as can be 
seen in Figure lb. Isospin violating processes are small but visible, such as the puj mixing 
seen in the pion formfactor. The situation for the meson is more involved. The OZI rule 
would prohibit decays into pions. However, a violation of this rule by about five percent 
is observed in the three-pion channel. (Some experimental data near the resonance 
have been left out in Fig. lb to emphasize the tu-contribution given by the solid line). 
The (f) meson decays mainly into OZI allowed channels such as K^K~ or KlKs (see Fig. 
Ic), but these channels also couple to the p and u mesons. The measured cross section 
of the annihilation into kaons includes an interference between all three vector mesons. 
(For more details see ref. |]14|). Nevertheless the meson still dominates the data in the 
e^e~ — s> KK channels. 



2.2 Vector Meson Dominance update 

In the region ^/s < 1 GeV vector meson dominance (VMD) describes the experimental 
data of e+e" annihilation into hadronic channels very well. We briefiy summarize here the 
results of ref. [|14| based on an effective Lagrangian which combines chiral dynamics and 
vector meson dominance. This approach starts out from a "bare", stable vector meson 

o ° o 

with a mass my slightly larger than the physical mass. Its propagator is Dv (? ) = 

(g^— my)~^. The coupling to the pseudoscalar mesons turns Dv into the physical vector 
meson propagator in vacuum, 

Dv{q^) = bv{q')+ bv{q')'nv{q')Dv{q') 

= ^ ' (9) 

g2_ my -ny(g2) 

with the complex vector meson self energy Uy- The decay width of a given vector meson 
is related to the imaginary part of Ily by 

myFy = — Imny(g^ = my), (10) 

o 2 

where my =my +Re ny(g^ = my) represents the squared physical vector meson mass. 



2.2.1 The p meson 

For example, the coupling of the p meson to the two-pion continuum leads to 



Imn^(g2 



9' 2 I, 



Am. 



2 \ 2 



q' 



0(g2-4m2), 



where g = 6.05 is the pTrvr coupling constant and = 139.6 MeV the charged pion mass. 
The real part of Up can be derived from a dispersion relation. It has the form 



Renp(g2) =cg2_ 



9' 



247r2 



q^ Q{q^,ml) - Am\ 



(12) 
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with the function Q{q'^, m?) given exphcitly in ref.|T^. The subtraction constant c = 0.16 



is determined by a fit to the isovector p-wave vrvr scattering phase shift, which also fixes 
the bare rho meson mass 0.81 GeV. Note that in this approach the physical rho 

mass in vacuum, rup = 0.77 MeV, is already shifted downward from the "bare" mass mv- 
One may think of this "bare" rho meson as being a quark-antiquark bound state before 
turning on its strong coupling to the vrvr continuum. ( In the l/N^ counting scheme, the 
bare rho meson mass is leading order in N^, whereas the coupling to the vrvr continuum 
via the complex self-energy Ily introduces l/N^ corrections.) 

2.2.2 The u meson 

The relation (p!0|) can be used to determine the imaginary part of the omega meson self 
energy by the decay into three pions (vr"'", vr", tt"), as follows: 

2 

Imn^(g2) = 1 1 dE+ dE.[plpl - (p„ ■ p^f] \Fs^\'. (13) 

Here E± and p± are the energies and momenta of the two charged pions in the final state, 
and the integration extends over the kinematically allowed region in the (£'+, E_) plane. 
The amphtude represents the sum of the direct u Sn decay and the Gell-Mann, 
Sharp, Wagner (GSW) process in which the omega meson converts first into vrp, followed 
by the decay of the p meson into two pions: 

i^3. = -TT + 2.4f E Dp{q - P.). (14) 

Jit In i 

Here we have introduced the pion decay constant fj^ = 92.4 MeV, and Dp is the rho 
meson propagator. The four-momentum pi refers to that of the first pion produced in 
the primary u —>■ np transition of the GSW process, and the index i denotes the three 
possible charges of that primary pion. The determination of the coupling constants in 
eq.(n) is described in detail in ref. [Q. The real part of H^^ can in principle be obtained 



by a dispersion relation. However, we choose to absorb this real part into the physical 
uj meson mass. This is justified since the momentum dependence of Re H^ does not 



influence the results in any significant way [0 , given the small omega meson width. 
2.2.3 The meson 

The self energy of the (j) nieson consists of three parts: 

= iifp^K+K- + n^^^o^o + H^-+3^ . (15) 

The last term violates the OZI rule, but despite the small (puj mixing angle 6 ~ 0.05 
it contributes about fifteen percent of the total decay width. It is given by eq.(|T3|) 
multiplied by the squared (pu mixing angle. The self energy parts coming from the decays 
into KK channels are of the same form as the vrvr self energy of the p-meson and connected 
with the latter through SU(3) symmetry: 



ReU^^K+K-{q^) = Ckq'^ - ^;^[q'^Q{q\mK+) - AmK+\ , (16) 
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Im n. 



Q 



1 - 



4m 



K+ 



(17) 



For n^^^o^o the same expressions hold with the charged kaon mass mK+ replaced by 
m^o. The subtraction constant ck = 0.11 has been determined in Ref.||14||. Comparing 
the theoretical widths with empirical ones we find that SU(3) symmetry works well and 



gives the same universal coupling constant g within ten percent WM for all couplings of 
the vector and pseudoscalar meson octets. 



2.2.4 Current Correlation functions and VMD 

We can now establish the relation between the reduced electromagnetic CC correlation 
function 11 and the vector meson self-energies Ily. Our starting point is again the effective 
Lagrangian of the coupled system of pseudoscalar mesons, vector mesons and photons as 
described in ref. [^. This approach includes not only a direct conversion of a photon 

o 

into a vector meson proportional to eg / Qy, but also a term which involves the photon 

o 

interaction with pseudoscalar meson currents. The couplings gy are related through SU(3) 

o o o _„ 

symmetry as Qp'.Quj'- 9^= 1:3:^. The same relation holds for gy, the strong couplings, 
where we identify gp = g. The result for n(g^), expressed in terms of Uylq"^), when 
iterating the vector meson self-energies to all orders is as follows: 




-ny(g2 



(1^ 



This can easily be verified by direct diagrammatic analysis (see Fig. 2) using the frame- 
work of ref. [n]. Here 7y(q'^) represents the photon - vector meson vertex with inclusion 
of vertex corrections due to pion or kaon loops: 



ay — 



g2 



(19) 



In case of the p meson we have ap = gj 1.1, the ratio of the pTrvr and p7 couplings. 
In the limiting case of "complete" VMD where all of the photon-hadron interaction is 
transmitted via vector mesons, we have ay = \. In the uj and channels the information 
about the direct 7 37r vertex is very limited and we therefore assume ai^ = a^ = 1. 
This gives an excellent description of the e^e~ —>■ tt+tt^tt^ cross section 
The imaginary part of 11 can now be easily evaluated: 



Imn(g^ 



E 

V 



lmlly{q'^ 



where 



Fv{q') 



[1 - ay)q'^- rky 



o 2 
m 



V 



-ny(g2 



(20) 



(21) 
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For example, in case of the rho meson, Fp{q'^) is simply identified with the pion form 
factor: 



9p; 



(22) 



where again g {dp) is the coupling of the p meson to the pion (photon) and Hp is the p 
meson self energy specified in eqs.(|Tl||l^). 



2.2.5 High energy limit 

In the high energy region the measured correlation function approaches the asymptotic 
plateau predicted by perturbative QCD: 

R{s) = -i^Imn(s) = V rfy fl + — ) e(s - sv) (23) 

S y \ IT / 

for s larger than a characteristic scale sy, where dp = 3/2, = 1/6 and d^ = 1/3. This 
behaviour is clearly seen in the p-meson channel where the low mass (two-pion) region 
of the spectral function is evidently given by p meson dominance in the form of eqs. (pOD 
and (P^D. For s > Sp ~ 1.5 GeV^ contributions from 4 and 6 pions take over and reach 
the asymptotic plateau above s ~ 2 GeV^. 



2.3 QCD sum rule analysis 

The basic idea of QCD sum rules is to match two representations of n(g^) in the region of 
large spacelike = —q^ ^ 1 GeV for each one of the channels V = p, uj, (j) (we drop the 
channel index and always refer to a specific channel in this section). One representation 
is in terms of a twice subtracted dispersion relation: 

n(g^) = n(0) + cq' + '-f ds ^.^"^^['^ .. , (24) 

TT J s^[s — — le) 

where the vanishing photon mass in the vacuum requires n(0) =0. The second represen- 
tation is in terms of the operator product expansion (OPE): 



127r 



n(g^ 



d 

TT 



(25) 



/iV 

(at a scale p which is commonly chosen around 1 GeV). Here the coefficients Ci^2,3 incor- 
porate the non-perturbative parts coming from the QCD condensates. In the isovector 
channel one finds 



Co 



rP 
^2 



—3(ml + 



TT 



1 + 
2 

3 vr 



-47r' 



(0|«s(^7m75A"m - J7M75A"rf)2|0) 



+ -{Q\as{u-ipru + d^p\'^d) gT'^A^glO) 

q=u,d,s 



(26) 
(27) 

(28) 
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Here 0:5 = 47r/ ( 6 In (Q"^ /-^qcd)) the running coupling constant (with Aqcd = 140 MeV 
and 6 = 11 — (2/3)A^/ = 9). We can set ci equal to zero because its contribution is 
proportional to the squared mass of the light quarks and therefore negligible. For the 
isoscalar channel there is only a sign change between the up and down quark condensates 
of the first term in coefficient C3. 

The situation is different in the meson channel because of the large strange quark 
mass nis. The coefficient cf is not negligible any longer and C2 is larger than 02'"^. One 



finds ^ 



1 + 



"1 



—6m, 



TT 
OS, 
TT 

{0\as{s-f,^,X''sf\0) 



q=u,d,s 



For the gluon condensate we use 



TT 



1.2- 10^2 GeV' 



(29) 
(30) 

(31) 
(32) 



the value determined from charmonium sum rules JT^ . The four-quark condensates which 
contribute to C3 are usually factorized by assuming that vacuum saturation holds, i.e. that 
they are proportional to (0|gg|0)^. Then the coefficient C3 is the same for all three channels. 
While this is a valid approximation in the large Nc limit |TB[, possible deviations from 
vacuum saturation must be taken into account. We introduce a parameter k > 1 for this 
purpose and write: 

Cs = -^KTT'as{0\qq\Of. (33) 

Here it is still assumed that the four-quark condensates with any two different flavours 
vanish. While this is not rigorously correct, the error can be estimated from pu mixing 
||17|| . Its smallness suggests that this assumption may be appropriate. With the commonly 
accepted values (0|mm|0) ~ {0\dd\0) ~ (0|ss|0) = (0|gg|0) ~ (-250 MeV)^ one needs 
K = 2.36 in order to obtain the same value of C3 as originally proposed by Shifman et al. 
p!7| . Note that Hatsuda et al. |jT9| use a larger value 



\qq\0) 



-280 MeV)^ in their 



approach. We also mention here that strongly increased values of C3 as suggested in refs. 
pO| , pT| would not be compatible with our analysis. 

The QCD sum rule approach compares the dispersion relation of eq.([2^) on the left 
side with its OPE representation (eq. (|25|) , right side). We use our calculated 



127r 
S9v 



lmUvis)\Fvis)\'Q{sv - s) + dy [l + — ] Qi-s - Sy 



as 



TT 



(34) 



to determine the left side. Since the OPE on the right side of eq. ( pSf ) is truncated at order 
we apply the usual Borel transformation [|T^, in order to improve convergence: 

d 



f{M') = ^ hm 

{Q2/n = Al^^const.) 



-too , n — 'oo , 



;n-i) 



iQ 



2\n 



dQ^ 



(35) 
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After Borelization the QCD sum rule for each R^^'^ becomes 



1 



roc 

1271^0(0) + / dsR{s)e' 
Jo 



Co + 



Cl 



+ 



C2 



+ 



C3 



2M 



(36) 



Given R^'^\s) of eq.(|3^) this represents a consistency condition which should hold inde- 
pendently of the Borel mass parameter Ai within the range of applicability of the QCD 
sum rules, i.e. for Ai larger than a certain minimal value. Note again in the vacuum 11(0) 
vanishes, but we have kept it here since it will not vanish in matter. 

In Fig. 3 we compare the "left side" of eg . (|36|) calculated from eg . (|3^) with the "right 
side" as a function of the Borel mass Ai . There is excellent agreement in all three channels 
between both sides of the sum rule, down to ^ 0.8 GeV where the convergence of the 
OPE series (|25|) starts to break down and higher order terms become important. We 
have checked that our condensates agree with the "optimal" values, using the fit program 
MINUIT. It determines the condensate parameters to be 



C3 



0.04 GeV^ 
-0.071 GeV^ 



(37) 
(38) 



together with Sp_^ = 1.56 GeV^ and 



C2 
C3 



0.07 GeV^ 
-0.10GeV^ 
-0.07 GeV^ 



(39) 
(40) 
(41) 



taking = 2.2 GeV^ 



We also confirm previous results which use the simplified ansatz for the spectra, 



as 



+ 1 + ^ e(s-sy 



(42) 



with T„ = ^Tuj 



2.4, viip^^ 



0.77 GeV, s 



1.5 GeV^ and Ta, = 0.79, 



1.02 GeV, 



S0 = 2.2 GeV . Including the vector meson decay widths does not make a substantial 
difference, at least not in vacuum: after taking the Borel transform their effects are 
suppressed as (myFy/A^^)^. This explains why the large difference of the p and u meson 
decay widths does not affect the OPE. However, for fine tuning purposes it is essential 
to use the complete spectrum (PW). For the in- medium sum rule the use of a realistic 
spectrum is mandatory as discussed in section |[ We also note that the approximate 
degeneracy of the p and u meson masses emerges naturally in the OPE. This is one of 
the traditional successes of QCD sum rules. 



3 Correlation functions in baryonic matter 

We now replace the vacuum |0) by a state \^{p)) of matter with given baryon density p. 
The CC correlation function to leading order in p has been written down before [^, . 
In this chapter we briefly summarize the steps which lead to this first order term, and 
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we also describe its iteration to all orders in density. The basic input into this scheme 
is the forward Compton amplitude for scattering of a timelike, virtual photon on the 
nucleon. Our task is to calculate this amplitude using an effective Lagrangian approach. 
The resulting spectral function in matter will then be further analyzed using a QCD sum 
rule with in-medium condensates. 



3.1 Basics 

In the present paper we restrict ourselves to the case of matter at temperature T = 0. The 
matter may be moving with a constant four- velocity v^. We choose a frame of reference 
in which matter as a whole is at rest by setting -u = 0. The four-momentum = (tu, q) 
that enters in the CC correlation function is from now on understood to be taken in this 
particular frame. The generalization of the CC correlation function to the case of baryonic 
matter is 

U,,{u;,q;p)=t J dt J d'x e'^'-'^^-^ {n{p) \rUt,x)jM\^ip)) ■ (43) 



For the vacuum case (p = 0) we have seen that at squared virtual photon masses q 



2 



UJ 



— g'^<lGeV^, the spectrum of 11 starts at = 4m^ in the isovector channel, at 
= 9m^ in the isoscalar channel and at g^ = 4mj^ in the (p channel (apart from small 
(f) ^ Stt contribution which also starts at g^ = 9m^). This spectrum displays the vector 
mesons as the dominant doorway states which decay into pseudoscalar mesons. In matter 
we have new processes involving nucleons, such as 

7* AT ^ ttN, ttttN, KA, iTS, etc. (44) 

which start at g^ > or at the threshold for kaon production. In addition there are 
nuclear particle-hole excitations at very low energies which need not be considered here. 

For a nuclear Fermi gas with equal number of protons and neutrons, Fermi momentum 
Pf and density p = 2p|,/(37r^), the density dependent matrix element {Q{p)\0\Q{p)) of 
any operator O is given by p3| : 



{n{p)\d\n{p)) = {o\d\o) + 4 / i%^An{p)\o\n{p)) + ... , (45) 

J\p\<PF {27ry E{p) 

where \N{p)) is a proton or neutron state (including spin) with momentum p and the 
normalization 

{Nip)\Nip')) = {2nf^6'ip-p'), (46) 

with E{p) = y/p'^ + M'^. The dots in eq. (^5D indicate higher orders in density. 
When applied to the CC correlation function in matter, eq.(^) becomes 

n^.(u;, q; p) = n;-(u;, q) + ^ I — ^— -T^.(u;, q;p) + ... , (47) 

J\p\<PF {■^TT) J^yP ) 
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where the first term UJ^ = Ili_iu{p = 0) is just the vacuum CC correlation tensor (|T|) 
discussed previously, and 



+00 



T,,{u;,q;p) = i dt x e'^'-'^^-^N [p) |Tj^(t, f)j,(0)| iV(p)) 



(48) 



is the forward Compton scattering tensor for a virtual photon of energy uj and momentum 
(f on a free nucleon. The tensor T^i, has the standard decomposition 



T, 



%<iy \rj. , ( p-q \ ( p-q \ T2 

- 9^i,^\Ti+ \p^, \ \p^ ^qu ' 



(49) 



-J y- ■ q^ ' / \ q^ ' / 
which defines the amplitudes Ti and T2. It is often more useful to introduce the amplitudes 



rj,(L) 



{p ■ qf 

p2 q2 



To 



(50) 
(51) 



and reexpress T^^, as 

T,, = -(Pl)m.^^''^ - (^t)m.T(^), (52) 

where the matrices Pl and Pt act as projectors on longitudinal and transverse parts with 
the properties: Pi = Pl, P^ = Pt and PlPt = Pt Pl = 0. For a nucleon at rest with 
= (M, p = 0) the tensors Pl and Pt are given by 



Pl 



Pt = - 







(53) 



3.2 Low density limit 

At low baryon density p the p dependence of the Compton tensor in eq.(P7|) can be 
neglected and we have 



n^,(^,g1p) ^ Iil^{q)+pT^,{u,q-p = Q) + 



(54) 



nr;(g) - P UPL),uT^''\u;,q,p = 0) + (Pt) ,.T^^\u; , q, p = 0) 



We introduce the notation T^^{u!,q) = T^^{u!,q;p = 0) etc. from here on. The reduced 
CC correlation function n(a;, q; p) = gfj,i,Ii^^ /?> analogous to eq.(^ becomes 



n(^,g;p) 



n^^^(g)-p(V)(.;,g) + ^T(^)(^,g: 
n™=(g) - p (t,{uj, q) + ^T2(^, q)\ , 



(55) 



to leading order in p. These expressions have a direct analogue when looking at the 
leading density dependence of the vector meson self-energy tensor: 



^v{^^q'^p) = ^T^''{q ) + pTvN{^,q 



(56) 
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where n™^'^'^ = ny^(p = 0) is the vacuum contribution and Ty'j^{uj,q) denotes the (off- 
shell) vector meson-nucleon amplitude analogous to the Compton tensor T'^'^ . With the 
help of the projectors Pt and Pl we can separate this self energy into two parts: 

u^y^iu, q; p) = [nr(g') - p4'i]{PL),u + [nr(g') - pTf;](PT),., (57) 

(L T) 

where Ily is the scalar self-energy function in the vacuum and Tyj^ are the longitudinal 
and transverse VN scattering amplitudes analogous to the Compton amplitudes T^-^'^^ of 
eqs.(0,0). 



3.3 Vector meson propagator and current correlation function 
in medium 

Using the projectors Pl and Pt it is now a straightforward exercise to decompose the 
tensor structure of the in-medium vector meson propagator. We find 



(PlT" 



q- 



VN 



+ 



[Pt) 



(T) 
VN 



(5^ 



The leading medium corrections, iterated to all orders, are therefore simply obtained 
by the replacement Ily^{q'^) — > Ily'^ — pTy'j^^ in the transverse and longitudinal parts, 
respectively. 

The structure seen in eq.(|58|) translates directly into the in-medium CC correlation 
function (^3]). Using diagrammatic techniques analogous to those in the vacuum one 
derives: 

(59) 



W^^icu, q; p) = n(^)(^, q; p)(P^)^'^ + n(^)(.;, q; p){Pj 



with the longitudinal and transverse parts 



n(^'^)(c^,g1p) 

2 

V 9v 



nr(r)-pr&^'(i^,« + 



w,q-p) 



V 



q2_ _n™^(g2) +pT^^P[u:, q ] 



(60) 



and 



'y'y'^\uj,q;p) = ay - 



W^^-pT^^p{uj,q) 



(61) 



This result generalizes eqs.(|l^,|T^). The imaginary parts of 11*^^^ and become: 



V 9v 



VN 



[1 - ay)q^- rhy 



q2_ _nvac + pT/^'^) 



■ VN 



(62) 



This is the in-medium generalization of eq. (|20|) , with 11™'^ replaced everywhere by IVy^ — 

P J-VN ■ 
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It is instructive to recover the low-density limit for which one finds: 



1 



n^.(a;, q- p) = W;:{q') - " Mq") T\,1j{Pl),u + T^'^(Pt);.. Fy{q') + ... (63) 

V 9v 

where Fy{q^) are the form factors given in eg. (PT|) . We can now identify the relationship 
between the Compton scattering amplitude and the VN amplitude by comparison with 
eq-O): 

^ ' ' ''^ '^^ >y (64) 



V 9v 



or 



(65) 



V 9v 

In the limiting case of a real photon with = u;^ — g ^ = we have Fv{q^ = 0) = 1 and 
find the simpler relation 



T^''^\oo=\q\)=Y. 



4'^P{u^=\q\) 



9l 



(66) 



Furthermore, 11™'^ (g^ = 0) = and ImT*^"^) = for real photons so that we obtain, to 
leading order in the density p: 



and therefore 



Imn^^.(u; = \q\\p) = -plraT^'^\uj = \q\) {PT)i,v, 
lmT^i{u = \q\) 



9'v 



UJ 

Aira 



(T^Ar(t^), 



(67) 



(6^ 



where a^N is the total photon-nucleon cross section averaged over proton and neutron 
and we have used the optical theorem 



eHmT^'^\uj = \q\)/uj. 



(69) 



3.4 Vector mesons at rest 

From now on we focus on the simpler case g = in which the vector mesons are at 
rest. Note that this is also the situation corresponding to back-to-back emission of an 
electron-positron pair. In this case we have 

T(^) = -^g^^T>^'^(u;,q = 0) = T,{u;,q = 0) = T^''\u;,q = 0) = T^''\u;,q = 0). (70) 

The reduced CC correlation function in the low density limit (eq.(^)) now becomes 

n(^, q = 0; p) = W^^iuj') -pT{uj) + ... , (71) 
and for the vector meson self energy we find 

Uv{uJ, q = 0; p) = nr (c^') - pTvn{uj) + ... , (72) 
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where Tt 



VN 



Ty^ denotes the (off-shell) vector meson-nucleon amplitude for 



(f = 0, analogous to the Compton amplitude Tioj"^). The full in-medium propagator 
for a vector meson at rest, simplifies to = Dy = Dy = and 



{my-%uj^) - pTvn) 



(73) 



The reduced CC correlation function in matter which follows from eq.(|59|) is now: 



V 9v 



VN 



o 2 



-Wy^-{tU^)+pTvN, 



with 



7y(u;;p) = av 



nr(^')-pTv^(c^) 



(74) 



(75) 



Its imaginary part enters directly in the back-to-back dilepton production rate and it is 
also the key input for the "left side" of the QCD sum rule: 



Imn(ci;, q = 0; p) 



E -rim 

V yy 



(1 — ay) tu"^— m 



V 



-ITy-^u^) + pTvNiuj) 



(76) 



We have thus established the relationship between the isovector CC correlation function 
in matter and the (complex) off-shell vector meson-nucleon amplitudes. 



4 Vector meson-nucleon scattering amplitudes 
4.1 Effective Lagrangian 

In order to determine the low energy VN scattering amplitudes, we use the hadronic 
effective Lagrangian of ref. [0] which combines chiral SU(3) dynamics with VMD. This 
approach has been successfully applied for vacuum processes. The extended Lagrangian 
incorporates nucleons, hyperons (and also A's in intermediate states ) with pseudovector 
meson-baryon interactions: 



Ub = FtT (57^75[m^ B]) +DtT (57,,75{m^ , (77) 

where B ($) represent the SU(3) matrix fields of the baryon (pseudoscalar meson) octet 
and 

u'^ = (d^^ - te[Q, $] A^) . (78) 
This effective Lagrangian has been tested successfully in the SU(3) coupled channels 



approach of refs. ||24|| . Here Q denotes the quark charge matrix, is the photon field 
and = 92.4 MeV is the pion decay constant. The second part of ■u'^ gives rise to 
the Kroll-Ruderman couplings of the photon with the interacting system of baryons and 
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pseudoscalar mesons. The parameters F ~ 0.51 and D ^ 0.75 satisfy the constraint 
qa = F + D = 1.26. The vector meson-baryon interactions are then introduced by the 
minimal couphng scheme which substitutes eQA^ by the vector field where 



^ pi' + uj'' ^ 

-p^ + cj^ 
\ V20^ J 



(79) 



In the KroU-Ruderman part this leads to 

U^B = (Ftr(57^75[[V^<l>],5]) + Dtr(57^75{[V^<l'],i?}+)) . (80) 



4/. 

The direct coupling of the photon to the baryons becomes: 

C^B = e (tr(57^[g, B]) - ii{B^^B)ii{Q)) (81) 



The second term vanishes, because the trace of the m— , d— and s— quark charge matrix 
is zero. When determining the direct vector meson-baryon coupling by the replacement 
eQA^ gV'^/2 this term survives and we find: 

C% = I (tr(57,[r^ B]) - tiiBj.BMVn) ■ (82) 

This form implies the SU(3) symmetry relation g = gp^ = \gujN = —^dcpN- We also 
include corrections due to anomalous VN tensor couplings with Hp = 6.0 and = 0.1: 

In the strict SU(3) limit we set = and recognize that there is no direct meson- 
nucleon coupling. The (pN interaction in this channel arises entirely from kaon loops. 

The interaction vertices obtained from the effective Lagrangian are summarized in 
Fig. 4. The minimal coupling scheme for the vector mesons makes sure that we correctly 
approach the limit u ^ keeping q = 0. 

Using this effective Lagrangian we now calculate the p, u and meson-nucleon scat- 
tering amplitudes Tvn{^), with the vector mesons at rest (g = 0). The imaginary parts 
of these scattering amplitudes are then evaluated using standard Cutkosky rules. The 
real parts are determined by a subtracted dispersion relation at g = 0: 

ReTyAr(t^) = /y H V / du —— -. (84) 

TT JO U^(n^ — UO'^j 



The subtraction constants ly are fixed by the Thomson limit at = 0. Using eq. (|66D 
we have Ip = 1^/9 = —g'^I^AMN) and = 0. Evidently, processes with low thresholds 
dominantly infiuence the low energy region of the real and imaginary parts of Tyjv. The 
convergence of the dispersion integral is slow, however, because ImTvN increases linearly 
with energy. Therefore its high energy part also infiuences ReTyTv- 
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Our hadronic effective Lagrangian is expected to be valid up to a momentum and 
energy scale of order 1 GeV and we expect formfactor suppressions at higher energies, in 
particular when the mesons are off-shell. Indeed, in the vacuum the VMD approach seizes 
to be valid at a scale s = sy, and perturbative QCD takes over for s > sy- We introduce 
formfactors for vector meson-baryon couplings: 

FvEik') = ^1^. (85) 

where k is the four momentum, my the mass of the vector meson, and Ay is the cutoff 
parameter. We choose Ay = 1.6 GeV, values close to those commonly used in Boson ex- 
change NN potentials. At the pseudovector pion-nucleon vertices we employ the empirical 
axial form factor 

with ~ 1 GeV and assume that a slightly larger A^'^ = 1.2 GeV applies for kaon 
couplings to baryons. We will discuss the sensitivity of the results to variations of the 
momentum scales A. 



4.2 The pN amplitude 

For the p meson the dominant contributions to the imaginary part of the scattering 
amplitude come from the inelastic processes pN nN, pN ttA ttttN and pN — ^ 
ujN — > titittN . The corresponding one-loop diagrams are illustrated in Fig. 5. The 
propagating baryon in the loops of this figure can either be a nucleon or a A isobar. In 
the limit of large baryon mass only the diagrams. Figs. 5(a,c,e,f,m), survive. First we 
omit the form factor (^^(A;^) and calculate the imaginary parts of the loops using the bare 
coupling qa. For the vrA^ loops we get 



g\n{LJ,0,m^) 



^7) 



where the function 7i is given (for uj > 0) by: 



H{uj, A, m) 



247r/2 



UJ 



Ay 



2cuA)^ 



+A^{lQuj^ - Sm^) - SA^o; + 4A' 



A(a;2 - 4m2)l 
'2u ro;2-4A2)2 



5 m 



4A^ 



If we include diagram 5 k, which does not contribute to the imaginary part of Tp^-, we can 
also evaluate the real part in the heavy baryon limit since all divergences cancel due to 
current conservation. We then find for the vrA^ loop contribution to Re Tp^ in this limit : 



ReT^ri^) 



2mi^ 



3m^ 



4m^ ^ 4m^ 



UJ' 



UJ^ 




ml — uj"^ 0(m7r — a;) 



m. 



^9) 
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This corresponds to the result obtained by using eq.(|87D as input in the dispersion relation 
(P^). Of course we have to take relativistic corrections and the axial form factor into 
account. Calculating the vrA^ loops fully relativistically and using Gyi(A;^) given in eq.(|86D 
instead of gA we end up with lmTl;]^\uj) as shown by the short-dashed line of Fig. 7a. 

Note that for the ttN loops the function ?i(ct;,0,m^) greatly simplifies. The situation 
is more involved when we consider the vrA loops. We find in the heavy baryon limit 

lmT^jf'\oo) = 2gln{uj, Ma - M^v, m^). (90) 

Here Ma = 1.232 GeV is the A-isobar mass. A fully relativistic calculation would require 
using the Rarita-Schwinger formalism which introduces ambiguities. Here we only take 
the axial form factor into account by replacing qa by a dipole ttNA transition formfactor 
G[^\k'^) with A^"* = A^. Adding this result to the vrA^ loops contribution gives rise to 
the long-dashed curve in Fig. 7a. 

For the high energy part of the spectrum we have found the box diagram 5(m) to be 
very important. Its fully relativistic evaluation gives 

lmT;,';,^\co)=g\liu,m^), (91) 

where 

^^^^^v) ^g^^4 M^(M^ + a;)(m2 -m2. + 2A;o^-a;2^2 ^^^f 



Here g^jfm = 1.2 and we use \ k\ = ^ X{M]\f + oo, mv, Mn) / {2{Mn+uj)) and ko = y my + k"^. 
The Kallen function A is defined by 

X{x, y, z) = {x^ -{y + zf) {x^ - {y - zf). (93) 

Including the axial formfactor and adding this to the contributions of the nN and vrA 
loops gives rise to the solid curve in Fig. 7a. 

One might expect that processes pN — > KS, KK also contribute to T^at. We give the 
result from KT, and KK loops in the heavy baryon limit, setting mx = 'mx+ = m^o: 



ImT;;^"'""^(^) = -iD- FY niu, - M^, tuk) + ttt(^ + 3F)^ Hiu, Ma - M^, ttik) 

o 24 

(94) 

Here M\ = 1.116 GeV and Ms = 1.191 GeV are the masses of A and S hyperons. Note 
that SU(3) factors reduce this amplitude already by about one order of magnitude with 
respect to the vrA^ loops, and a further reduction comes from the much smaller phase 
space and relativistic corrections. We can therefore safely neglect those contributions. 

The real part of the pN scattering amplitude is then determined by using ImTpAr (solid 
line of Fig. 7a) as input for the dispersion relation (0). The result is shown if Fig. 8a. 
The decomposition of ReTp^ into subprocesses is also given in this figure. 

4.3 The (J)N amplitude 

For the 0-meson the same contributions from the KA and KTi loops as for the p-meson 
appear, but now larger by a factor of two: 

ImT^Niu;) = ^{D-Ffn{u;,Mj:-MN,mK) + ^{D + 3Ffn{u;,Mj^-MN,mK) (95) 
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Note that from SU(3) symmetry the (pN amphtude is already four times smaller than 
the pN amplitude with inclusion of ttN loops only. A further reduction comes from the 
smaller phase space because of the higher KN threshold. Therefore the imaginary part 
and consequently also the real part (shown in Fig. 8c) of the (f)N amplitude is much 
reduced as compared to the pN case. This is evident by comparing Figs. 7a and c. The 
smallness of ImT^^^v implies that effects of regularization by meson-baryon formf actors 
have little impact on the real part, at least in the region around the physical (f) mass. We 
show results in Fig. 8c with inclusion of a kaon-baryon dipole formfactor using a cutoff 
A^'^"* = 1.2 GeV. We mention however that diagrams such as Fig 5(1, m), with K* and K 
mesons exchanged in the box are not taken into account in the present calculation. Very 
little is known about such processes, but they may not be negligible. 



4.4 The ujN amplitude 

The KroU-Ruderman term and the coupling to two pseudoscalar mesons do not exist for 
the oj meson. Consequently no such diagrams as in Figs. 5(a,b,c,e,f,h,i) contribute to the 
imaginary part of the ujN amplitude. The diagrams. Figs. 5(d,g,j) vanish in the heavy 
baryon limit. They contribute very little to the scattering amplitude even with the large 
coupling = 17. The dominant contributions come from processes where the uj turns 
into a Tip system which then interacts with the nucleon. The corresponding box diagrams. 
Fig. 5(m) (with the role of p and uj interchanged) give: 

ImTi]}(u;) = Sfif^ J(u;,mp), (96) 



where X(a;,mv) is defined in eq. (p2D. At low energy the box diagram Fig. 5(1) becomes 
important since the rho mesons in the box interact strongly with the nucleons due to the 
large anomalous tensor coupling with Kp = 6. We find for this term: 

, ^(2). . 9lp.9\^ + ^P? \l\ (^' - mlf {AMI -< + 4Mjva; + u^) 
Vml^^{u)- ^^o_xo 777 — --:T7Z3 ZTo .9,0, .7 • 1^0 



128 7r/2 (M^ + uY {ml - ml - tu^ + 2ujk 



using |/| = ^\{Mpf + cu, m^r, M^) / {2{Mn + t^)) and Iq = y "ml + where A is the Kallen 
function given in eq. (p3|) . 

Including also the contributions from vertex corrections in which the u couples on one 
side directly to the nucleon and on the other side to the vrp system, we show the resulting 
imaginary part of the uN amplitude in Fig. 7b replacing by the axial form factor and 
multiplying Kp by the monopole form factor (^) . At low energies ImT^^^v is small but rises 
strongly at higher energies. The high energy behaviour is sensitive to the meson-baryon 
form factors. This unavoidable model dependence also translates into ReT^p^ (see Fig. 
8b) when using the calculated ImT^^ in the dispersion relation (Q). We should also point 
out, however, that the behaviour of ImT^M can at least partly be checked against data. 
For example, the mechanisms that determine the long-dashed curve in Fig. 7b also enter 
in the s-wave part of the irN — > uN cross section. We find indeed satisfactory results for 
this cross section close to threshold. 
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4.5 Effective vector meson - nucleon scattering length 

In order to gain further insights into the basic features of the VN amphtudes taken at 
g = 0, it is instructive to introduce complex effective scattering lengths, 

= A in^^A -JvN{i^ = my), (98) 

4:71 [Mn + my) 

defined for "on-shell" vector mesons at rest. These scattering length parameters should 
of course be interpreted with care since the vector mesons themselves are unstable. Nev- 
ertheless the real and imaginary parts of avN give a direct impression of the strong 
interaction shifts and widths that the vector mesons would experience in nuclear matter 
at low density. 

The following effective scattering lengths are obtained for our standard scenario, with 
input specified in section 4: 

ttpN = (0.04 + U.62)fm, (99) 
a^TV = (3.34 + i 2.1) fm, (100) 
a^jY = (-0.01 +z 0.08) fm. (101) 

We observe that the pN scattering length is strongly dominated by its imaginary part 
which describes the inelastic pN — > tcN, tttcN etc. channels, whereas Re is persistently 
small, more than an order of magnitude smaller than Im Cp^r. In essence, this OpTv causes 
a strong inelastic broadening, but almost no mass shift in matter. This is at striking 
variance with almost all previous analyses which have simply ignored the large imaginary 
part of ttpN- 

The situation is quite different for the ujN scattering length. There one finds a large 
(attractive) real part of a^^N which is mainly driven by the processes. Fig. 5(l,m). Thus 
the u meson mass is indeed expected to be shifted downward in matter. At the same 
time one observes a substantial inelastic broadening. 

Finally, the comparatively weak (j)N interaction is reflected in the small a^at, with 
almost vanishing real part and a moderate Ima^jv which comes primarily from the inelastic 
channels (pN ^ KA, KE. 



The results, eqs.( p9| - p!01| ), are of course subject to some model dependence related 
to the form factors at the meson-baryon vertices that enter into loop diagrams. Fig. 5. 
Changes by 20% of the cutoffs A in these form factors induce variations at the 30% level 
in the real and imaginary parts of avN- For example, increasing the cutoff A^"* of the 
dipole ttA^A transition form factor from 1 GeV to 1.3 GeV will change the pN scattering 
length to ttpN = (—0.02 + il.4)fm. It is thus uncertain whether the very small p meson 
mass shift in matter is upward or downward, but the primary feature of a large Im apN 
remains in all cases. We have also examined corrections due to vrvr rescattering in the 
basic pN -k-kN box diagram (see Fig. 6) and found them to be small. Such mechanisms 
produce attraction, but their effects on the p mass shift in matter are still weak (less than 
5% at p = po) 
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5 Spectral densities and QCD sum rules in medium 



5.1 In-medium spectral functions 

Combining our results for the VN scattering amplitudes with the relations developed in 
chapter 3 we are now able to compute the in- medium correlation function (^,0). We 
insert Tvn into eq.(|7^ and display the spectral function in the form — (127r/a;^) Imn(u;, p) 
for convenience, the one corresponding to eq.(||). This is done separately in each of the 
vector meson channels for different baryon densities p. 



5.1.1 Isovector spectral function 



In the isovector channel (see Fig. 9a) we observe the following characteristic features 
as compared to the vacuum: 

- The imaginary part of the in-medium rho meson self energy is strongly increased by 
the large Im Tp^r. This leads to a huge width in the CC correlation function, with the 
consequence that the rho meson resonance broadens and flattens. 

- The spectrums starts now already at u = m^r, the threshold energy at which single 
pion production becomes possible. This generates a strong continuum in the correlation 
function below the p resonance. We recall that the vacuum spectrum starts at a threshold 
LJ = 2m^ coming from p —> tt+tt" decay. 

- The shift of the p meson pole position in matter is only marginal, given the small 
real part of the pN scattering length. The in-medium rho meson mass is determined by 

^lip) =K [^TH) - pTpNimD] . (102) 
With our standard scenario we find 

m%P) = mliO) (^1 - 0.005-^j , (103) 

a very small shift at nuclear matter density, po = 0.17fm~^. 

It is instructive to study the real part of the in-medium meson propagator, 

D^{uj,q = 0;p) = -. (104) 

u;^-m^-(U;^{u^)-pTMu^)) 

The zero of KeDp determines the in-medium mass. The result is shown in Fig. 10a. One 
clearly observes that, with rising density, the rho meson has little chance to survive as 
a quasi-particle. The large inelastic width in matter diffuses the well developed quasi- 
particle pole structure seen in the vacuum. The actual downward shift of strength in the 
spectrum in matter is caused entirely by the strong increase of the imaginary part of the 
in-medium self energy. 



5.1.2 Isoscalar spectral function 

For the u meson the situation is different (see Fig. 9b). It starts out as a narrow resonance 
in vacuum. Both real and imaginary parts of T^^n are large. The u meson mass shifts 
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downward and its width increases by an order of magnitude at p = po as compared to the 
free width, but it survives as a quasi-particle in matter, at least up to densities p ~ po- 
This is clearly seen in the real part of the in- medium uj meson propagator (see Fig. 10b), 
the zero of which moves downward with increasing density while the pronounced resonance 
pattern still persists. 



5.1.3 The spectral function 

Because of the relatively weak ipN interaction, the meson peak moves very little (see 
Fig. 9c). At p = Po its width increases by about one order of magnitude over its very small 
decay width in free space due to the opening of the inelastic channels (f)N — > i^A, KTj. 
The spectrum now start at the threshold for those processes, but the well developed (f) 
resonance remains clearly visible at p = po. 



5.2 QCD sum rules in baryonic matter 

Our calculated in-medium CC correlation function at g = is now used as input for the 
dispersion relation which represents the left side of the QCD sum rule at finite baryon 
density. We then examine the consistency with the operator product expansion, the right 
hand side of this sum rule. 

As in the vacuum case described in section 0, it is useful to work with the quantities 



/?(^)(c.;p) 



127r 
9v 



Imlivi^iQ = 0,p) 6(sy 



+ d\ 



1 + ^ ) Q{u^ 



Sv) 



(105) 



in each of the three vector meson channels {V = p, u, 0). The calculated spectrum Imlly 
(see eq.([76D) is valid up to some (density dependent) scale sv{p) which still needs to be 
fixed. At higher energies beyond sy the perturbative QCD result takes over; here we 
assume that the coefficients dy do not change with density. After Borel transformation 
one obtains the in-medium QCD sum rules in a form analogous to eq.( 



9v 



rfa;2i?W(cu,p) e--'/-^= 



Co + 



Cl , C2(p) , C3(p) 



+ 



2M^ 



(106) 



The first term on the left hand side corresponds to the Landau damping, the vector meson 
analogue of the Thomson scattering limit, with Tvn{0) = gy/^^N for V = p, uj. For the 
meson this term vanishes. 

The coefficients ci, C2 and 63 include the density dependence of condensates which 
have already appeared in the vacuum sum rule, as well as new condensates which only 
exist at finite baryon density. For the isovector and isoscalar channels we have B: 



~P,Ll) 



^2 



-0 ) 



Cl 



0, 

2 ; 



(p) = cr(p) + 2vr^Ar'^M;vp. 
10 



5r(p) = cr(p) - yvr^^r'^M^p, 



(107) 
(108) 

(109) 
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where Cg'"^ are the same as in the vacuum. For we use 0: 



o vr 



^{oi^g^-'g.M-^MS^p, (110) 



where m'^^ = 750 MeV is chosen for the nucleon mass in the chiral hmit [^] . The 
contribution from the chiral condensate, {Q\muUU + rridddlQ) = {OlmuUU + miidd\0) + pa ^ 
with the nucleon sigma term ~ 45 MeV, is small and can be neglected. 

The in-medium four-quark condensate determines 03''^ (p). Its density dependence is 
not very well under control and requires some discussion. The basic question is to what 
extent ground state saturation, {fl\{qqy\Q) ~ {Q\qq\Q)'^ , is approximately vahd in matter. 
This factorization of {{qq^) holds in the mean field approximation where one neglects any 
excitations that can be reached by acting with qq on the ground state. In a nuclear medium 
as well as in the vacuum, scalar multi-pion excitations of the ground state are expected 
to be relevant, however. By analogy with eqs. (pSj) and (Rl) we parameterize: 



1 fi 

Mp)m,7^^'qrMp)) = -{^{pmi.^'qmp)) = -fn{pmp)\qqMp))\ 

where the density dependent factor /t(p) now represents the in-medium deviation from 
simple factorization. Ground state saturation means k = 1. We now assume that the 
leading density dependence comes from the in- medium condensates {Q{p)\qq\fl{p)) = 
{0\qq\0) + aN p/imu + rrid), while the low-energy part of the spectrum of scalar excitations 
that enters in k{p) > 1 is approximately independent of p (see also refs.[|, ^). We 
therefore choose /t(p) = /t ~ 2.36 as in the vacuum and find 



cr = vr^ as UO\qq\Or - —^{0^10) (HI) 

to leading order in density p. This may be a crude approximation which overestimates the 
density dependence. At least the tendency for {{qqY) to decrease with increasing density 
is without question [^], since the difference between vector and axial vector current 
correlation functions, ~ (IjaJaDI, is proportional to {qq)"^ and tends to zero at 

the point of chiral restoration. 

As pointed out in ref. 0, further in-medium contributions to eqs.( |107] - |109|) arise from 
matrix elements of the form {N\q'^^^D^^...Dfj,^j^-^q\N) with gluon fields entering in the 
gauge covariant derivatives D^. These matrix elements are proportional to the moments 

Al = 2 dx x'^Mx) + q{x)]. (112) 

JO 

of quark and antiquark distributions which can be determined, at some renormalization 
scale of order 1 GeV, from deep-inelastic lepton scattering on a nucleon. We use values 
as in ref.[|] and set A\'^'^ = 0.9, A^'^'^ = 0.12. Additional corrections of higher twist 
operators are smaller than the twist 2 contributions fS^, ^ and discussed in ref. . 



22 



For the (j) meson the coefficient C2 is dominated by the strange quark condensate. Its 
in-medium change is given to leading order in density by 

{n{p)\rssMp)) = {0\ss\0) + y (113) 

with y = {N\ss\N)/{N\uu\N) ~ 0.2. We then find 

ct = ci ~ct = ct (114) 

ctip) = ct{p)+A7c^AlMNP, (115) 
on 

ctip) = c^(p) - yvrM^^^p, (116) 

where C2(p) and c^ip) are obtained from 02,3 of eqs.( ^0| - |3l| ) when replacing the vacuum 
|0) by \Q{p)). We use the values = 0.05 and = 0.002 of the strange quark moments 
as in ref. 0. The resulting density dependent changes of the right hand (OPE) side of the 
QCD sum rule for the meson with respect to the vacuum are small. 



In Figs. 11 we compare the "left hand side" with the "right hand side" of eq.( 106 ) 
for the p, u and channels at normal nuclear matter density p = Po = 0.17fm~^, as a 
function of the Borel mass A4. For the energy scales sy in eq.( |105|) which separate the 
resonant parts of the spectrum from the high-energy QCD continuum, we use 

Sp = 1.56 GeV^(l -0.1 p/po), (117) 
= 1.65 GeV^(l -0.3 p/po), (118) 
= 2.2 GeV^(l- 0.01 p/po), (119) 

in order to match the region of asymptotic Borel masses. 

The overall consistency between calculated in-medium spectral functions and the den- 
sity dependent operator product expansion is evidently quite remarkable. We should 
compare Figs, lla-c with the corresponding vacuum sum rule analysis, Figs. 3a-c. Note 
that the difference between the p = po and p = cases is large for the p and u meson 
channels, whereas it is marginal for the (p channel. 

5.3 Discussion 

5.3.1 Comparison with previous work 

Hatsuda and Lee P| have inserted the simplified parameterization 



i?(^)(c., p) = J^;.S{cu' - m*y'{p)) + dy{l + ^ ) e{u;' - 4(p)) (120) 



in the "left hand side" of their in-medium QCD sum rule and fitted J^y, rriy and Sy. 
Apart from their choice of k = 1, their "right hand side" is the same as eq.( |106|) . Their 
optimal fits gave 



m 



Pi ^ 



1 - (0.16 ±0.6)-^, (121) 



rrip^^ Po 
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= 1 - (0.15 ±0.05) — , (122) 
Po 

= 1 - (0.24 ±0.07) — . (123) 

J^p,uj Po 

Here J-'p = 127r^m^/(7^ is the vacuum pole strength for the p meson, and an analogous 
expression holds for the u meson. For the in-medium meson mass they find 

^ = 1 - (0.15 ±0.5)?/-^, (124) 

Po 

with y = 0.1 — 0.2. Jin and Leinweber find relations consistent with eqs.( |121[ - |124| ). If 
we use K = 2 the results ( [L21[ - |124D change moderately within the errors indicated. 



The fact that the in-medium masses m* ^ both appear to decrease by (10-20 %) at 
p = Po in comparison with the vacuum masses rrip^ ^ has been interpreted as an indication 
of BR scahng [|l|. In contrast, we find in our explicit calculation that the low-energy 
interaction of p and u mesons with surrounding matter differ strongly. The naive ansatz 
( |120|) remains reasonably justified in the isoscalar channel. There we have shown that, 
while the effective width of the u meson at p = po increases by about an order of magnitude 
over its (small) free width, the in-medium u meson can still be identified as a quasi-particle 
with a well-defined effective mass that roughly resembles the fit, eq.( p.21| ). This is not the 
case, however, for the p meson. The inelastic collision broadening in the isovector channel 
is so strong that it completely dominates the low-mass spectrum, whereas the pole mass 
of the rho meson at p = po changes very little. The average spectral mass rh in this 
channel, defined by 

drops roughly as m* in the parameterization ( p.20|) , but this can obviously not be inter- 
preted as a decreasing p mass in matter, as we have pointed out by analyzing the zero of 
the real part of the rho meson propagator in Fig. 10a. 

Asakawa and Ko [0 have used a more realistic p meson spectral function than the 



schematic zero width form ( |120| ) and repeated the QCD Sum rule analysis with inclusion 
of a density dependent rho meson width. They have argued that a mass shift consistent 
with BR scaling still remains to be added. Our results for the p meson do not support 
this conclusion. 

The strong broadening of the p meson spectral function in matter is also found by 



Chanfray et al. I^Gf. They include corrections of higher order in density induced by 
intermediate propagating pions. Our present calculations have focused on a consistent 
calculation of vector meson self-energies to leading order in density. We find that uncer- 
tainties in the (small) p meson mass shift related to vertex form factors and high energy 
contributions to the spectral functions are in fact larger than the effects of higher order 
corrections in pion propagators. 
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5.3.2 Effective scattering amplitudes 



The main difference between the present approach and previous work is our explicit cal- 
culation of the complex vector meson- nucleon amplitudes at g = 0. In the low density 
limit the question about possible shifts of vector meson masses reduces to a statement 
about the corresponding effective pN , ujN and scattering lengths. In previous discus- 
sions the large imaginary parts of these scattering lengths have simply been ignored. For 
example, Koike and Hayashigaki |22| deduce a real pN scattering length apN — 0.45 fm, 
so that the relation m*^ — rn?p = — 47rp(l + nip/ Mi^)apN implies a downward mass shift of 
about — 45MeV at p = po- Our calculated effective scattering length (|98D is qualitatively 
different. We find a large imaginary part hxvapM ~ 1 fm which reflects the strong influence 
of the pN — i> nN, vrA inelastic channels. The real part Re OpTv, on the other hand, is about 
one order of magnitude smaller than ImOpTv- It may have different signs depending on 
details of cutoff procedures, but this is almost irrelevant in view of the large magnitude 
of the imaginary part. 



Hatsuda, Lee and Shiomi [|T9[ have pointed out that the low-density limit 5m 



—pTvN is expected to break down already at a fraction of nuclear matter density po. This 
is indeed the case, given that the magnitude of Tvn is large for V = p, uj. It is necessary 
to use the full vector meson propagator (^), with the self-energy Ily = —pTvN iterated 
to all orders, and this has in fact been done when deriving the spectral function (|76D. 
Additional higher order corrections in Uy as calculated e.g. in ref. do not lead to 
substantial changes for p<po- 



6 Summary and concluding remarks 

(1) Guided by the chiral meson-baryon effective Lagrangian with inclusion of vector 
mesons and electromagnetic interactions, we have calculated the spectral distributions 
of current-current correlation functions in each of the p-, u-, and 0-meson channels, both 
in the vacuum and in baryonic matter. We observe perfect consistency with QCD sum 
rules in all cases. 

(2) Our results do not suggest a universal scaling law for in-medium vector meson masses. 
We find almost no mass shifts for p and (/)-mesons in matter, whereas the u meson expe- 
riences a substantial (attractive) shift. 

(3) The pN and ujN effective scattering amplitudes are governed by large inelasticities. 
For the p meson channel the inelastic broadening moves a large amount of the spectral 
strength down to the continuum below the free p meson resonance. This broadening is so 
strong that the p meson is unlikely to survive as a quasiparticle at densities much above 
Po, the density of normal nuclear matter. In contrast, the u meson, with its small width 
in vacuum to start with, can still be identified as a quasiparticle in matter at p = po. The 
(pN interaction is altogether weak. Consequently, the meson is expected to have almost 
no mass shift. Although its width at p = po increases by about an order of magnitude 
over its free width, it still persists as a relatively narrow resonance. 
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(4) To the extent that the changes of vector meson spectra in matter are driven by an 
effective chiral Lagrangian, they represent the dynamics and symmetry breaking pattern 
characteristic of low energy QCD. These features enter at the same time in the operator 
product expansion part of the QCD sum rules. Whether they directly reflect a tendency 
toward chiral symmetry restoration at high density is less obvious. 

(5) Our calculations indicate that the lifetime of the ui meson in matter at p = po decreases 
to about ~ 1.5 fm/c. The production of "slow" uj mesons in a nuclear environment 
can therefore make in-medium changes of the uj meson spectrum visible. Such processes 
can be explored with the upcoming HADES lepton pair spectrometer at GSI. A reaction 
of particular interest is Tx~p un in heavy nuclei with subsequent decay u e~^e~ and 
a kinematic cut focusing on low u meson momenta ||39|, E0[ . 



(6) For the lifetime of the meson we still expect ~ 10 fm/c at p = po which ex- 
ceeds nuclear dimensions. We mention that once the effective in-medium and K~ 



masses are introduced according to ref. pT[, a slight downward shift of the meson is 



expected (by 1-2 % at p = po), and its hfetime decreases to about 6 fm/c. 
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Figure Captions: 



Figure 1: The ratio a{e^e^ hadrons)/(j(e~'"e~ — >• fJ^^fi") as a function of the total 
cm. energy. The isovector, isoscalar and the meson channels, corresponding to the 
currents (3-5), are shown separately in Figs, la, b and c, respectively. The data are taken 
from refs. [30-35]. In the channel the high energy data sample summarizes KK + nir 
final states. The solid lines are the results of the VMD model. The dashed lines represent 
the perturbative QCD limit. 

Figure 2: Diagrammatic representation of the full photon propagator (wavy line) and 
of the full vector meson (V) propagator. 

Figure 3: Comparison of the Borel transformed right side (operator product expan- 
sion, OPE) and left side (dispersion relation, DR) of the QCD sum rules for the isovector 
(1=1), isoscalar (1=0) and meson channels, as a function of the Borel mass Ai (see 
eq.(|36D). The solid line is the result using our spectral density which describes the ex- 
perimental data. The short dashed line is the simplified ansatz using a delta function for 
each vector meson resonance. 



Figure 4: Feynman rules for the vertex functions derived from the effective Lagrangian 
of the coupled system of vector mesons, pseudoscalar mesons and baryons. Examples are 
shown here for interactions involving p mesons and pions (including the npuj coupling). 
Solid lines: nucleons; double solid lines: A-isobars; dashed lines: pions; vector mesons are 
drawn as "zig zag" lines. The generalization to flavour SU(3) is straightforward. Here S 
and T are the A — ^ spin and isospin transition operators and is the isospin operator 



for the A (see ref. 42 



Figure 5: Processes contributing to the pN scattering amplitude. The dashed lines 
represents pions, the wavy line describes the p meson. The incoming baryon (solid line) 
is always a nucleon while the intermediate baryon can be either a nucleon or a delta. 



Figure 6: Two-loop contribution to the pN scattering amplitude including tttt inter- 
actions. 



Figure 7: Imaginary parts of the VN scattering amplitudes in the isovector (a), 
isoscalar (b) and meson (c) channels. Solid lines: total results. Dashed lines show 
decomposition into selected subchannels. 



Figure 8: Real and imaginary parts of VN scattering amplitudes in different chan- 
nels. Solid lines: imaginary parts used as input for dispersion relations to generate the 
real parts shown as long-dashed lines. Fig. 8a also shows separate contribution from 
pN ttN (with qa = 1) and pN nN + ttA (with full axial formfactor). 

Figure 9: Spectra of the current- current correlation functions in different channels, 
shown at zero (dashed lines), half (long dashed) and normal nuclear density (solid lines). 
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Figure 10: Real part of in-medium vector meson propagator. The lines correspond to 
zero (dashed), half (long dashed) and normal nuclear density (solid). 

Figure 11: The Borel transformed "left side" (l.s.) (from dispersion relation, dot 
dashed) and "right side" (r.s.) (from OPE, solid line) of the QCD sum rule (see eq.( |106| )) 
at normal nuclear density. For comparison we also show the vacuum results of Fig. 3. 
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